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Abstract 

Spinning black three-branes in type IIB supergravity are thermodynamically stable up 
to a critical value of the angular momentum density. Inside the region of thermody- 
namic stability, the free energy from supergravity is roughly reproduced by a naive 
model based on free M = 4 super- Yang-Mills theory on the world-volume. The field 
theory model correctly predicts a limit on angular momentum density, but near this 
limit it does not reproduce the critical exponents one can compute from supergrav- 
ity. Analogies with Bose condensation and modified matrix models are discussed, and 
a mean field theory improvement of the naive model is suggested which corrects the 
critical exponents. 
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1 Introduction 



Low-energy physics on the world- volume of N coincident D3-branes [1] placed in flat 
ten-dimensional spacetime is dictated by M = 4 super- Yang-Mills theory (SYM) with 
gauge group U (N) [2] . The gauge coupling and the closed string coupling are related 
by 9ym = 27r 9s- 

A cluster of coincident D3-branes can also be regarded as a black three-brane of 
ten-dimensional type IIB supergravity (SUGRA) [3]. The geometry has characteristic 
curvature scale R = (inggNa' 2 ) 1 ^. Since the validity of supergravity depends on 
having curvatures small on both the string and Planck scales, both N and gy M iV must 
be large. Large 't Hooft coupling, A = gy M iV ^> 1, is the opposite limit to the one 
treatable in perturbation theory, so if supergravity does capture certain aspects of 
gauge theory, it is of strongly coupled gauge theory. 

One of the earliest attempts to exploit the relation between supergravity and gauge 
theory as alternate low-energy descriptions of D3-branes was the computation of the 
entropy [4, 5]. The results disagree only by a numerical factor: in the microcanonical 
ensemble (which seems the most natural from the point of view of gravity), 

ssugra = 
ssym = (4/3) 1/4 ssugra ■ 

In (1), s and e are entropy and energy per unit world- volume. In supergravity, s is 
computed using the Bekenstein-Hawking formula, and e is computed as the ADM mass 
density above extremality. In super- Yang-Mills theory, s and e are computed using the 
free theory (gyM — 0). Conformal invariance dictates the e 3 / 4 behavior but not the 
coefficient. The disagreement in the coefficient is sometimes quoted as 4/3 rather than 
(4/3) 1 / 4 ; this is the result of working in the canonical rather than the microcanonical 
ensemble. Given that the gravity and gauge theory calculations are valid in the limits 
of large and small 't Hooft coupling, respectively, the puzzle is not so much that there 
is a numerical discrepancy in the coefficient but rather that the discrepancy is finite 
(that is, it does not depend parametrically on N or the 't Hooft coupling A) . It is as if 
the number of degrees of freedom in supergravity were exactly 3/4 the number in the 
free theory. Indeed, it was pointed out in [4] that the discrepancy could be cured by 
excluding the gauge bosons and one species of gauginos from the free field counting. 
Because the gauginos transform in an irreducible representation of the R-symmetry 
group, this prescription breaks rotational invariance. A more likely scenario (recently 
challenged in [6]) is that the coefficient on the entropy interpolates smoothly between 
finite weak and strong coupling limits. It was shown in [7] that at least the leading a' 
correction to supergravity perturbed the coefficient toward its free field theory value. 
In section 2 we will extend the supergravity calculation to the case of the black three- 
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brane with angular momentum in a plane perpendicular to the brane. The geometry 
was constructed explicitly in [8], following the work of [9]. This spinning three-brane 
has the same quantum numbers as N coincident D3-branes with a density of -R-charge 
on the world- volume equal to the angular momentum density j (the /^-symmetry group 
is precisely the rotation group SO (6) in the dimensions transverse to the brane). Fol- 
lowing the spirit of [4], we suggest that the supergravity geometry actually teaches us 
the strong-coupling thermodynamics of large- iV super- Yang-Mills theory The salient 
feature of that thermodynamics is an instability for j > N 2 T 3 (up to numerical factors 
of order unity) . 

It seems impossible to recover this prediction by a standard perturbative gauge 
theory calculation. Zero coupling is in a sense the worst place to be expanding around, 
because in the presence of charged massless bosons there cannot be any "voltage" for 
the charge density j. In section 4 we will suggest a formal approach for regulating the 
free field partition function in the presence of such a voltage. We do not claim that 
our approach invalidates the free field theory conclusion that massless bosons cannot 
support a chemical potential; however, it may be a better starting point from which to 
understand the strong coupling regime. In particular, we shall see that near j = the 
regulated free energy reproduces the supergravity results for thermodynamic quantities 
like specific heat and isothermal capacitance, up to factors of order unity. Furthermore, 
it also predicts thermodynamic instability for j > N 2 T 3 , again in agreement with 
supergravity up to a finite factor. The critical exponents are predicted incorrectly by 
the regulated field theory model, but we suggest that this might be fixed by a mean 
field treatment of interactions. 

Our attention will be focused, both on the supergravity side and in the regulated 
field theory model, on the phase whose region of thermodynamic stability includes the 
states of small T and zero j. We believe that the boundary of stability is actually a 
phase boundary beyond which another stable phase emerges. Our speculation is that 
in this phase the D3-branes spontaneously separate. The motivation for this guess will 
be explained at the end of the paper, but the bulk of it is devoted to three goals: 1) to 
demonstrate the existence of a boundary to the region of stability of the small j phase; 

2) to examine the critical properties of the thermodynamics near the boundary; and 

3) to see how the thermodynamics might be reproduced in field theory. 

2 Spinning three-brane thermodynamics 

Since the rotation group 5*0(6) transverse to a three-brane in ten dimensions has 
rank three, it is possible to give a three-brane three independent (commuting) angular 
momenta. We will use the solution constructed in [8], which has only one angular 
momentum, but there should be no difficulty in principle using the results of [9] to 
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extend to the general case. The metric is 

,9 1 / , ,9 n.\dr 2 2£rnCOsho; . ,, 

ds 2 = —= i-hdt 2 + dx 2 ) + Jf °—— r — sin 2 6dtd(p 

Vf y ' V 1 h r 4 A/ Y ( 2 ) 

+ r 2 {Ad6 2 + A sin 2 #d0 2 + cos 2 9dn\ 

where 



/ = 


1 + 


r 4 sinh 2 a 


r 4 A 


A = 


1 + 


e cos 2 e 


A = 


1 + 


e ! r^ 2 sm 2 


h = 


1 - 


4 

r 4 A 


h = 


-A 


r 2 r 4 1 



(3) 



and the dilaton is constant, as for the spherically symmetric case. Fixing the number 
iV of units of five-form flux leads to the constraint 

Nk 

Tq sinh a cosh a = R 4 = - 5 ^ 2 (4) 

where 

2k 2 = 16nG = (2n) 7 g 2 a' 4 (5) 

is the gravitational constant of type IIB supergravity. 

We define the energy of this solution as the ADM mass minus the mass of a super- 
symmetric, non-rotating three-brane with the same five-form flux: if V is the spatial 
world- volume, then 

= E = M ADM - M 
' V V (6) 



— — (j csch 2a + coth 2a - 1 J 



K 

where we have used the BPS formula for the tension of the extremal three-brane, 
Mq/V = Ny/n/n. The angular momentum conjugate to can be read off from the 
asymptotics of [10]: 

jxy ^3 

j = — = —£r 4 cosh a . (7) 

V Kj 

Note that e has dimensions of (length)" 4 while j has dimensions of (length) -3 . Finally, 
the entropy 
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Ah/AG 2k 

s = — - — = —r H r cosh a (8) 

V K 
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is computed from the area A H of the event horizon, whose location r H is the unique 
real positive root of h = 0: r\ = \(sjl* + 4r 4 ~ ^) [Sl- 
it is possible in principle to eliminate £, r , and a from (4), (6), (7), and (8), 
and express s = s(e,j;N,R), where we regard (e,j) as variables parameterizing the 
phase space and N and R = {An g S N a' 2 ) 1 1 14 as parameters specifying the theory. The 
nontrivial dependence on the dimensionful parameter R makes it hopeless to try to 
reproduce this entropy function from Af = A SYM, because M = A SYM is conformal: 
it has no scale. However, it was observed for j = already in [4] that the -R-dependence 
drops out in the near-extremal limit (Madm — M )/M < 1 (a startling consequence 
is that all dependence on the 't Hooft parameter drops out as well, in the leading 
a' approximation to type IIB string theory which type IIB supergravity represents). 
For j = one is left with an entropy function s = s(e; N) which is unique up to an 
iV-dependent normalization given conformal invariance of the underlying theory and 
extensivity of S and E. For j ^ the entropy function again becomes independent of 
R in the near-extremal limit r ,£ <^ R, but it involves a function of the dimensionless 
ratio j 4 /e 3 that cannot be determined by scaling arguments. We will now determine 
it from supergravity. 

The near-extremal limit has a — > oo, and a is conveniently eliminated in favor of r 
using an approximation to (4), e 2a = ™* r i ■ One immediately learns e = §^2?"o> an d 



with a bit more work j = t\j -JiiNe/'SK and 

' 4 3- 3 / 4 V^ 

/ v /r + ^+v / X JJ X (1 + X) + X 



s = s(e,r, N) = 25/ ; 3 ' 3/4 ^ 3/4 = , ^ (9) 



where 
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The denominators represent the nontrivial new information we are gleaning from the 
spinning black hole solution. 

As a safeguard against algebraic errors, and to check the normalization of (7), we 
have computed the angular velocity Q at the horizon, fl is defined by requiring the 
Killing vector d/dt + Qd/d(j) to be null at the horizon. One finds 

The first expression for Q is exact; the second represents the near-extremal limit. In 
this limit, one can check from (9) that 

n = -W§Y- < 12 > 

(os/oe)j 

as required by the first law of black hole thermodynamics: dE = TdS + QdJ. 
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3 The limits of thermodynamic stability 



Our primary interest in this section is the thermodynamic stability of near-extremal 
spinning D3-branes. As a warmup, let us study the stability of non-rotating black 3- 
branes which are not required to be near-extremal. The entropy function s = s(e; N, R) 
has a rather complicated form when written out explicitly. It is simpler to refrain from 
eliminating the variable a from (6) and (8), since then we have 

e = — — — ( | csch la + coth 2a — 1 ) 

271 R (13) 
1 sinh 1 / 4 ^ [l6) 



T = 



ttR cosh 3/4 a 



In this simple example, thermodynamic stability is equivalent to positivity of the spe- 
cific heat at constant volume, Cy = (dE/dT)y. In manipulations of thermodynamic 
quantities, we will always hold V fixed, and we will usually drop the subscript V that 
indicates this explicitly. Also we will work with intensive quantities only. Inspection 
of (13) reveals that while e is strictly decreasing with a, T has a quadratic maximum 
at a = a c where coth 2 o; c = 3. This indicates a singularity, cy ~ 1 / \JT — T c , in the 
specific heat, since 

de de/da 

Cv = df = WJda- ' (14) 
and T ~ T c — [a — a c ) 2 near the critical point. For a < a c , cy < 0, so this region 
is unstable. Note that the scale of the phase transition is set by R, so it can't have 
anything to do with renormalizable M — 4 super- Yang-Mills theory. If it has any 
world-volume interpretation, it seems it would have to be in terms of an effective, 
possibly non-local theory like the Born-Infeld generalization of electromagnetism. On 
the supergravity side, the phase transition is much less arcane: it is the standard 
transition from positive to negative specific heat that occurs as mass is added to almost 
any charged black hole solution. The thermodynamic instability is real in the sense 
that Schwarzschild black holes really do get colder when you increase their mass. 

What we will find for the near-extremal limit of spinning D3-branes is a boundary of 
thermodynamic stability with the same \/y/T-T c singularity in a specific heat, but 
with the critical temperature T c set by the angular momentum density rather than a 
scale built into the theory. In other words, it is an effect which falls inside the putative 
regime of validity of the renormalizable super- Yang-Mills theory description of the 
world-volume. Spinning three-branes with T < T c are unstable, and the question we 
begin to address next is how one might see this boundary of stability in field theory. 

It is slightly more subtle to evaluate thermodynamic stability when the entropy 
function S = S(xi, . . . ,x n ) depends non-trivially on several extensive thermodynamic 
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variables X{. The criterion is sub-additivity: 



S{{\ Xl + (1 - \)x t }) > \S({ Xl }) + (1 - X)S({x i }) 



(15) 



for A between and 1. When S({xi}) is a smooth function, one can simply demand 
that the Hessian matrix [d 2 S/dxidxj] has no positive eigenvalues. 

Applying the above considerations to the entropy function (9), we see that the con- 
dition for thermodynamic stability is \ < 1/3- Equivalently, we may demand that the 
energy function 




(16) 



^2/3(2^)2/3 

derived from (9) has a positive definite Hessian. This also leads to % < 1/3, or more 
usefully 

(17) 



j 1 
- < -7=- 



Since the Hessian matrix [d 2 e/ dxidxj] can also be written as [dyi/dxj] where y,i = 
de/dxi, a positive definite Hessian means that the functions yi = Vi(xi, . . . ,x n ) are 
invertible — ie that we can Legendre transform freely among the various ensembles. 

The ensemble which will be most directly comparable to the naive field theory model 
developed in section 4 is the one where (s,j) are both traded for their dual variables 
(T, Q). We will call this the grand canonical ensemble because it is like having indefinite 
particle number and a chemical potential. The free energy in this ensemble is usually 
called the grand potential, and we will denote it as S. It is related to E by 



E = E + TS + m . 



(18) 



Written as a function S(T, Q, V) of the volume plus the intensive quantities temperature 
and voltage, the grand potential is forced by extensivity to be linear in V . Indeed, 
— £ = — S/V = P is just the pressure. Conformal invariance implies £ = —\e. In 
deriving the explicit expression 
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f(T,S2) = -^-JV 2 T 4 



(i + ^S^-v^ff) 2 )' 

(i + ^S(i-v /rr ?f)7 



(19) 



7T 



= —N 2 T A 



in 2 in 4 

1 + t?T 2 + 2t^T 1 



+ ... 



the following "equation of state" is useful: 



f 



1 +2n 2 j 2 /s 2 



(20) 
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The region of allowed Vl/T is restricted: 



which as the reader can check expresses the maximum value the right hand side of (20) 
can attain as a function of j / s. The boundary of stability (17) occurs at the maximum 
value of VL/T, so the (21) is not a good description of the boundary of stability: the 
equality is obtained at the boundary, but on both sides Q/T decreases.* It is possible 
to quote the stability condition in the canonical ensemble, where T and j are fixed: 
the condition is 

j = —^-=N 2 T 3 . (22) 

There is really only one critical exponent in the description of the thermodynamics 
near the boundary of stability. It characterizes the approach of £ to its value at 
criticality: at fixed temperature, 

Z = Z c + j c (n c -n)-A(n c -n) 2 -< + o[(n c -n) 2 } . (23) 

where A is some constant and, for the case at hand, 7 = 1/2. All other critical 
exponents follow from this one using scaling properties. In particular, the specific heat 
cq at constant voltage goes as (T — T c )~ 7 . With 7 = 1/2, this is the same 1 / y/T — T c 
behavior we saw in the non-rotating example at the beginning of this section. 



4 A naive field theory model 

Let us first recall the spectrum of M = 4 super- Yang-Mills theory. Besides the gauge 
bosons, which have no charge under the R symmetry group SO (6) ~ SU(4), there is 
a (color) adjoint of Weyl fermions A a in the 4 of SU(4) and an adjoint of scalars X 1 
in the 6 of SO (6). Under the subgroup SO (2) of SO (6) generated by <9/<90, half of 
the fermions (that is, AN 2 degrees of freedom) have charge 1/2 and the other half have 
charge —1/2, whereas of the scalars, two thirds (AN 2 degrees of freedom) are neutral, 
one sixth have charge 1, and the last sixth have charge —1 (note that we have included 
particles and anti-particles explicitly in our counting). All these particles are massless 
since we have not Higgsed the theory. 

The presence of charged massless bosons makes it impossible to apply a "voltage" Q 
of the SO (2) subgroup to the free theory, simply because one of the thermal occupation 
factors n±(p) = l/(e^l^ ±n ^ — 1) will become negative for \p\ < \Q\. The interactions, 
treated perturbatively in the 't Hooft coupling A, result in a thermal mass m ~ \f~\T 

"Thanks to F. Larsen and P. Kraus for setting mc straight on this point. 
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to lowest order. Thus it seems the critical value of the potential f2 is fi c = \f\T up to a 
factor of order unity: parametrically different from the supergravity result. Factors of 
A will also enter into thermodynamic quantities like isothermal capacitance (dj/dQ)x 
measured at j — = 0. There is nothing contradictory in supposing that in fact 
fi c = /(A)T for some function /(A) which approaches a constant for A — > oo and \/X 
for A — > 0; and likewise for the various measurable quantities like (dj/dfl) T Q _ Q - 

Positing an interpolating function for every quantity that disagrees between free field 
theory and supergravity may be in some sense the correct thing to do, but it is highly 
non-predictive, particularly since both the perturbation expansion in finite temperature 
field theory and the a' expansion in string theory are difficult to compute past the first 
nontrivial term. In hopes of establishing a better field theory understanding of the 
supergravity results, we will suggest in section 4.1 a method of regularizing the free 
field theory partition function so that it becomes well-defined at finite voltage. In 
section 4.2 we will see that this naive field theory treatment has some degree of success 
in reproducing the phenomena observed in supergravity. In section 4.3 we will try to 
improve the naive model using mean field theory. 



4.1 A regulation scheme 

In this section we will suggest a regulation of integrals over n±{p) which makes pre- 
dictions for the critical voltage Vt c and for (dj /dVi) T which agree up to factors of 
order unity with the supergravity results. It is not our aim to revise free field theory; 
rather, we hope that our procedure captures in some sense a better approximation to 
the dynamics of the effective excitations of interacting super- Yang-Mills theory. The 
factor of A^ 2 on the grand potential £ in (19) is an indication that there are N 2 degrees 
of freedom excited at finite temperature, even in the limit of strong coupling. We 
regard this as a clue that perhaps some adaptation of the free theory can be used to 
approximate the statistical mechanics of the interacting one. 

The divergent integrals that occur in free field theory with a chemical potential for 
massless bosons are of the form 

/„ *^3T = r (")? i V" = r(n)Li " (e " ) (24) 



f 

JO 



dpp"- 1 log(l - e""*) = -r(n) Li n+1 (e") . (25) 

where Li n is the polylogarithm function, which is defined by the sum in (24) when 
its argument is in the unit disk. This function can be analytically continued to the 
complex //-plane, and it is analytic except for a branch cut along the positive real axis. 
The imaginary part changes sign across this cut while the real part is continuous. The 
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principle value prescription for the integrals in (24) and (25) for real /x is to average 
the values for \i + ie and /i — ie. This amounts to replacing Li„ by its real part 3?Li n . 
One might be inclined to regard the principle value prescription as a species of infrared 
regulation, similar in spirit to differential regulation in its use of analytic continuation 
to define divergent integrals. But it differs from the use of an explicit mass as a 
regulator in that it yields a cutoff-independent free energy which is stable in a finite 
region including zero voltage. We will show that if both regularizations are applied at 
once, then by making the explicit mass small one smoothly recovers the massless case. 

Using this prescription on (24) one can calculate the charge density j for a given 
potential Q by integrating the thermal occupation factors. We will instead use (25) 
to compute the grand potential density £ = S/V = £(T, Q;N), from which all other 
thermodynamic quantities follow by differentiation. For a collection of several species 
of free massless charged particles in a potential fl, 



$ = T J ^ £ s t log(l - Si e-«ia-* n >) = ~T,*i LU(*e*» n ) . (26) 



The sum over i is over particle species. For each species, q is the charge and s = 1 
or —1 according as the particle is a boson or a fermion. For fermions, the principle 
value prescription is superfluous because the polylogarithm is evaluated at negative 
argument and so is already real. For free Af — 4 U (N) SYM, (26) specializes to 



N 2 T 4 



■3? 



7T^ 



-—N 2 T 4 



6Li 4 (l) + Li 4 (e-^ 2 ) + Li 4 (eH - 4 Li 4 (-e^ sz/2 ) - 4Li 4 ( 

4 2 Q 2 _ 1 tt A ~ 
3 + ^ 2 T 2 ~4^T* 



where we have used the identity (valid for positive integers n) 

i n ~ 2j ± 



[n/2\ 2/-(2i) 
Li n (en + (-1)" Li^e"") = E J^jjy''' 



ITT 



(n-1)! 



n-1 



(27) 



(28) 



The explicit sign on the imaginary term in (28) is the same as the sign of the imaginary 
part of ji. From (27) one can extract the charge density: 



n i ft 3 

T ~ 4^2 T3 



(29) 



As an aside, it is interesting to note that the two terms in (27) pertaining to the 
charged scalars individually have singularities in their second derivatives at the origin, 
but their sum is C°° . If we were to change the charge of the scalar or give it a mass 
without making the corresponding change to its anti-particle required by CPT, this 
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singularity would reappear. In a non-relativistic system of massive particles (that is, 
the usual Bose-Einstein setup), 



where we have again applied the principle value prescription to the integral for jj, > 0. 
dj/d/j, is singular at the origin and negative for fi > 0, indicating thermodynamic in- 
stability. Thus in the case of non-relativistic bosons, one recovers through the principle 
value prescription the standard conclusion that the disordered phase in a free Bose gas 
is stable only up to a charge density j = C(3/2)(mT/27r) 3 / 2 . 
Let us now consider the case of a relativistic massive boson: 



where we have included the anti-particle in the second term. We will regard the mass 
as an explicit infrared regulator of the massless theory. The integral is not one of 
the standard special functions. However, its singularity structure as a function of the 
complex parameter <£> is easy to determine: it is analytic except for branch cuts on the 
real axis for |<E»| > m. The branch cuts are discontinuities of the imaginary part. If 
one applies the principle value prescription for real $, one gets a continuous function 
j($) which is smooth except for a "kink" at 1$! = m, where the function itself remains 
finite but its first derivative is singular. In standard free field theory, this kink is the 
signal for Bose condensation. But if we continue past this kink using the principle value 
prescription, and if m is not too large, we find a new region where dj/dfi > 0. The 
size of the kink scales as (m/T) 3 / 2 . In the m — > limit, it disappears altogether and 
we recover a cubic curve qualitatively identical to (29) (it differs in exact coefficients 
because (29) includes the effects of fermions). The convergence of j(Q;m) to j($;0) 
as m — > is uniform. If the principle value prescription is used uniformly, physical 
quantities defined in a theory with an explicit mass as the infrared cutoff have a smooth 
behavior as that infrared cutoff is removed. Note however that the cubic curve (29) 
shows no sign of Bose condensation at $ — > 0. Its thermodynamic stability will be the 
subject of section 4.2. 




(30) 




(31) 
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4.2 Thermodynamics of the naive field theory model 



The reader who is distrustful of the regulation methods used in the section 4.1 may at 
least be willing to regard 



7T 



e = —n 2 t 4 



4 2 ft 2 
3 + ^T2 



1 ft 4 
I^T 4 



(32) 



as a motivated guess for the field theory partition function. But is it a good guess? In 
this section we will show that (32) exhibits a boundary of stability at a value of 
differing from the supergravity result (22) only by a finite factor, but that the critical 
exponent 7 is different. In section 4.3 we will argue that the correct critical exponent 
can be obtained through a mean field theory improvement of (32). 

The agreement up to factors of order unity of the first two terms in square brackets 
in (32) with the corresponding terms in (19) implies a similar agreement in the small ft 
limit of all thermodynamic quantities expressible in terms of second derivatives of the 
various free energies: such quantities can be related to entries in the Hessian matrix of 
£(T,ft), which is 



d 2 c 

dT 2 
dfldT 




- m 2 7r 2 T 2 



fiVVT 2 



1 

Q 



3n 2 T 

4 
3 

20. 
3tt 2 T 



n 

6^2 ■ 

2Q 
3tt 2 T 

3^2 . 



(l + 0[ft 2 /T 2 ]) forSUGRA 
( 1 + O [ft 2 /T 2 ] ) for field theory. 



(33) 

The upper left entry reflects the well-known 4/3 disagreement of the non-rotating case. 
The ft 4 /T 4 term in square brackets in (19) and (27) are of opposite sign, and so also 
some of the 0[ft 2 /T 2 ] terms which we have not indicated explicitly in (33) disagree in 
sign between supergravity and field theory. This is in fact the first symptom of the 
disagreement in critical exponents at the boundary of stability. 

It is already a nontrivial fact, however, that the regulated field theory predicts a 
boundary of stability at finite ft/T. Let us see how it arises. Thermodynamic stability 
requires that the Hessian matrix of £(T, ft) be negative definite. This is so when 



ft 2 

- < 271-4/-= - 1 

T ~ V VS 



(34) 



Note that this is a smaller value than the one we might guess by analogy with Bose 
condensation from the behavior of j(T, ft): the maximum of j(T, ft) for positive ft and 
fixed T is attained at ft/T = 2n/y/3, which is greater than the right hand side of (34). 
It is more enlightening to examine the "equation of state" 



J_ 

s 



3 ftl 



1 n 2 



4tt 2 T 1 + -2-0^ 

4vr2 T 2 



(35) 
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Figure 1: cjjasa function of Q for fixed T, in supergravity (solid line) and in the naive 
field theory model (dashed line). The horizontal and vertical axes have been scaled 
differently for the two curves so that the boundary of stability occurs at the same point 



Q c , and so that the two curves meet at Q — 0. 



Here the maximum of the right hand side does occur precisely at the boundary of 
stability 

By expanding £ in (27) at fixed T around the critical value of Q, one finds the 
behavior (23) with 7 = 0. No derivatives of £ become singular: for example, at fixed 
temperature, cq ~ (Q c — f2)~ 7 = const. As one can see from figure 1, cq smoothly 
approaches a maximum as Q — > Q c . 



4.3 Interactions and modified critical behavior 

The comparison between supergravity and the naive field theory model is a little like 
the comparison between superfluidity in 4 He and Bose condensation in an ideal gas 
[11, 12]. The free field treatment of Bose condensation predicts the location of the He 
phase transition to within 30%. But whereas in the free field treatment the specific 
heat has a finite peak at the transition temperature, in real 4 He it exhibits a log(T — T c ) 
divergence. The interactions modify the critical behavior. 

There is also an example of logarithmic corrections to critical scaling in an exactly 
soluble model of two-dimensional quantum gravity: the modified matrix models [13, 
14, 15, 16, 17]. In matrix models one typically has an expansion 

T ~ (analytic in A) + A 2 ~ 7 + (higher order in A) . (36) 

for the free energy. Here A is the cosmological constant, which is taken to zero to 
obtain large random surfaces, and 7 is termed the string susceptibility. The standard 
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c = 1 model has A 2 / log A as the leading nonanalytic term. One can introduce a 
(tr$ 2 ) 2 term, which can be regarded as an interaction among otherwise free fermions 
in a potential well. The positions of the fermions are the eigenvalues of the matrix $. 
If the strength of the interaction is tuned appropriately, it alters the critical behavior 
of c = 1 to A 2 log A (and, incidentally, changes the location of the critical point in the 
coupling space of the matrix model by a finite factor). This is like going from 7 barely 
negative to 7 barely positive. For c < 1, the shift is finite, and again changes the sign 
of 7: upon adding the interaction, 7 — > 7 = 7/(7 — 1). 

It is our hope that by including an effective description of the interactions in our field 
theory model, we may see a similar modification in the critical behavior — optimistically, 
one which will predict 7 = 1/2. 

Let us try to evaluate whether this is a realistic hope by considering a mean field 
theory treatment. One can suppose that the presence of a density j of i?-charge sets up 
a "self-field" contribution to the voltage. ^ In order to work in terms of dimensionless 
quantities (as conformal invariance demands), we assume a relation 

"tot " . /■ / j 



+ (37) 



T T 

where f(p) is some function to be determined. We obtain a consistency condition by 
replacing j in (37) with the integral over thermal occupation factors: 



(2tt) 3 Y e/3(|p1 ~^ tot) - S 

where as before % indexes particles, qi is the charge, and Sj = 1 for bosons and —1 for 
fermions. The integral is regulated using the principle value prescription. (38) boils 
down to replacing Q by Q tot in (29) and then plugging the resulting expression for j 
back into (37). Note however that the physical charge density must be obtained by 
differentiating the grand potential £ by the applied voltage Q. The grand potential is 
obtained by replacing Q by f2 tot in (27). Using the consistency condition one can write 
£ wholly in terms of Q and T. 

To shorten the formulas, let us introduce some notation for dimensionless quanti- 
ties: 

Jpf^ = P f= u — = ^tot. (39) 
Then the grand potential can be written as 

£ 7T 2 1 

- g(u) = g(uj tot ) = — + W ot - -30;^ , (40) 



jy2^4 a\ / i»V tot/ g 1 4— tot 32tt 2 



^This assumption may seem strange in view of the fact that SO(6) i?-symmetry is not gauged in 
the world-volume theory. It is however gauged in the supergravity. If the dynamics of supergravity 
are thought to be contained in the strongly coupled gauge theory, then it seems quite natural for two 
i?-charges to feel each other's influence. 
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and the consistency condition is 

u = uj tot - f(g'(u to t)) = wtot - / (k^tot ~ g^2 w tot) • ( 41 ) 

Primes will always denote the differentiation of a function with respect to its argument. 

The quantity p = du to t/duj can be regarded in analogy with OB/dH in magnetic 
systems as a permeability: p > 1 is the analog of paramagnetism, where the self-field 
contribution reinforces the applied voltage, and p < 1 is the analog of diamagnetism. 
Let po be the permeability at zero voltage. Modulo a constraint on the behavior of 
f(p) for finite p, we claim that 7 = 1/2 arises precisely when p = 2/3, and otherwise 
7 = or 3/2 according as po > 2/3 or p < 2/3. The constraint is that f(p) must be 
such that 

±g(u)gT(u) - 3~g'(cu) 2 > (42) 

for uj less than the value where u tot = uj t ot,c = 27T. This is a special point because 
g'(2n) = 0. We will refer to (42) as the convexity constraint because it is precisely the 
requirement that the Hessian matrix of £(T, fl) be negative definite. 

First note by explicit differentiation of (41) that p — 2/3 is equivalent to f'(0) — — 1. 
Let us assume po = 2/3. Then dui/duitot = at u tot = u tottC . Near u tot = cj to t,c we 
have 



uj c -uj = a 2 (cJ totiC - u tot ) 2 + a 3 (cJ totiC - cJ tot ) 3 + 0[(uJ tot)C - u tot ) 4 } 



(43) 



ff(Wtot) = ^0 + fl^tot.c - ^tot) 2 + S^tot.c - ^tot) 3 + 0[(wtot,c - ^tot) 4 ] • 

Plugging the first equation into the second, one obtains 

g(u) =g + gi(u c -u) + g 3 /2(uj c - ujf /2 + 0[(u c - uj) 2 } (44) 

which indicates that indeed 7 = 1/2. The reader interested in the details will have no 
trouble verifying 



UJ C 


= ^tot,c = 27T 






a 2 


3 

47T 


a 3 = 






_ ^ 




i 

2 


9o 


3 


92 = 




2vr 2 




2tt 


9o 


3 


9i = 





8tt 2 

1 (45) 
^ 3 = 47 

^3/ 2 = ^/|(3 + 2 7 r 2 r(0)). 

In (45) we have used the fact that f(—p) = —f(p) by symmetry, so all even derivatives 
of f(p) vanish at the origin. The critical charge density is given by 

= Pc = ~9i = y ■ ( 46 ) 
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For p < 2/3, du/dutot < at u tot = cj tot c . This makes it impossible to invert (41) 
to obtain u tot as a single-valued function of ui past the maximum value for u>, which is 
obtained at some value u to t = k>tot,m < ^tot,c- Because g'(u)tot,m) ^ 0, this scenario leads 
to 7 = 3/2. For p Q > 2/3, u) totiC could still be the critical point if (42) is not violated 
earlier, but there would be a linear term in u to t,c — k>tot in the expression for u c — u in 
(43), so g{oj) would be analytic in uj near oj c : 7 = 0, as we saw with no interactions 
included at all. So only at the special value p = 2/3 can we see the 7 = 1/2 behavior 
predicted by supergravity. The precise shape of g(cu) depends on the behavior of f(p) 
for finite p. But the location oj c = 2n of the boundary of stability and the prediction 
7=1/2 rely only on /'(0) = — 1 — again, modulo the convexity requirement (42). It 
may be entertaining for the reader to note that with po = 2/3 the discrepancy between 
the isothermal capacitance in field theory and supergravity is a factor of 8/9, which is 
numerically an improvement over the factor of 2 that we saw for p — 1 in (33). 

At this point we must ask if it is indeed possible to satisfy (42) for < u tot < 2n, 
because if (42) is violated for u to t < c^tot.c — 2vr, then the system becomes thermody- 
namically unstable "by accident" before one can come to the special point u to t = c^tot.c, 
and all our analysis will have been for nothing. As long as u is a strictly increasing 
function of u tot , (42) is equivalent to requiring that 



so a necessary condition (and also a sufficient one if du/du tot > can be shown) is to 
require (47) for < uo tot < 2ir. We do not have complete results on what restrictions 
this implies on the behavior of /(p). However, it is straightforward to verify that linear 
functions f(p) = 2 ii ^-p satisfy the requirement if and only if p — 2/3. For p < 2/3, 
linear f(p) results in 7 = 3/2 in the way described in the previous paragraph. For 
Po > 2/3, in fact the critical point occurs for c<j tot < u; to t,c — 27r, and 7 = by an 
analysis qualitatively identical to the case with no interactions at all. 

It is conceivable that through some peculiar singular behavior of f(p) at finite p one 
could arrange 7 = 1/2 and a smaller critical value for u tot . Of course, such a setup 
would seem extremely ad hoc unless one had a justification for the particular form of 
f(p) from an analysis of the interactions. Such an analysis is highly desirable at this 
point, at least for small p. It would be particularly amusing if the strong coupling limit 
A — > 00 managed somehow to tune p to 2/3. Indirect evidence for whether this is 
happening might be extracted from a study of a' corrections in the supergravity, along 
the lines of [7]. 




(47) 
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5 Conclusions 



We have investigated the thermodynamics of coincident D3-branes spinning in a plane 
orthogonal to the branes. Using the standard techniques of black hole thermodynam- 
ics we obtained an expression for the entropy which indicates thermodynamic stability 
only up to a special value of the dimensionless ratio of the angular momentum density 
and the temperature. We can write all thermodynamic functions in terms of inten- 
sive quantities because conformal invariance forbids the presence of any dimensionful 
parameters: any function of dimensionful variables can be re-expressed in terms of 
intensive quantities by scaling each variable by the appropriate power of the volume. 
The absence of a scale means that physical relations like an equation of state or the 
location of a boundary of stability must be expressible in terms of dimensionless ratios 
of intensive quantities. It also implies that only one critical exponent, 7, describes the 
approach of the world- volume theory to the boundary of stability. It for instance gives 
the scaling law ~ (T — T c )~ 7 at fixed voltage (as explained above, voltage Q is the 
thermodynamic variable dual to i?-charge). We have found 7 = 1/2 from supergravity. 

In the superconformal field theory which characterizes low-energy fluctuations of the 
branes, angular momentum is interpreted as electric charge under a subgroup of the 
/^-symmetry group. We have suggested an alteration of free field theory in which one 
can calculate very simply the critical value of the angular momentum density where the 
phase transition takes place. The calculation employs the thermodynamic ensemble 
where voltage rather than _R-charge density is held fixed. This ensemble is formally 
ill-defined when there are free massless charged bosons. The alteration of free field 
theory which we have suggested relies on a regulation procedure employing analytic 
continuation to render all the relevant integrals finite. The reader is entitled to question 
the validity of our regulation scheme: negative thermal occupation number is a serious 
pathology. Nevertheless, if supergravity is supposed to be teaching us about super- 
Yang-Mills theory at strong coupling, then spinning black holes show that the theory 
is stable up to voltages on the order of the temperature. A description of the interacting 
theory has to start somewhere, and our starting point, the regulated free field theory, 
enjoys a measure of success. It leads to an expression for the grand potential which 
exhibits a boundary of stability similar to the one observed in supergravity, and at a 
similar value of the charge density at fixed temperature. 

The regulated field theory wrongly predicts 7 = 0. In analogy to Bose condensa- 
tion and the modified matrix models, we have suggested that including an effective 
description of the interactions may change the scaling behavior. We have substan- 
tiated this claim by showing in mean field theory that if the three-branes respond 
"diamagnetically" to an applied voltage (that is, if the total voltage on the branes is 
reduced from the applied voltage by a self-field contribution) then the scaling behavior 
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can be modified to 7 = 1/2. One obtains 7 = 1/2 at a special value, /x = 2/3, of 
the "permeability." This value separates a region of 7 = behavior (/j,q > 2/3) from 
a region of 7 = 3/2 behavior (/j,q < 2/3). The regulated free field theory and the 
mean field treatment of interactions are not motivated by any deep understanding of 
strong coupling gauge dynamics, but rather by a desire to see in a simple way how and 
whether finite-temperature field theory can reproduce the thermodynamics of spinning 
near-extremal D3-branes. 

The analysis of this paper leaves us at a point similar to the state of matrix models 
of two-dimensional quantum gravity: we have a description of critical behavior near a 
boundary, but on the other side of the boundary our description breaks down. In the 
matrix models, the intuition was that eigenvalues were leaking out of the local potential 
well and falling down a V(X) ~ —A 4 hill. A clear description of the physics seemed 
difficult because it demanded a "non-perturbative" definition of two-dimensional quan- 
tum gravity. In Bose condensation a new stable phase (the Bose condensate) emerges 
on the other side of the boundary. The obvious analog of Bose condensation on D3- 
branes is Higgsing: some of the scalars develop expectation values. This suggests that 
the cluster of coincident branes should be splitting up, perhaps into a ring (condensed 
state) around a central core (normal state). A supergravity solution of this type may be 
difficult to write down explicitly. If it could be shown that there is a correlation length 
in the scalar fields on the uncondensed side of the phase transition which diverges at 
the boundary, we would regard it as evidence in favor of the Higgsing scenario. We 
hope to report on these issues in the future. 

Note Added 

After this work was completed, we received the paper [18] which has a tangential 
relation to section 2 of this paper. In [18] the Euclidean version of the spinning three- 
brane solution was considered. The thermodynamics is qualitatively different because 
the angular momentum and velocity are imaginary. The entropy function that follows 
from (3.24) of [18] can be obtained formally from our results by replacing ~^ ~\fx 
in the first expression for s in (9). 
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